Quantitative Atomic Spectroscopy for Primary Thermometry 
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Quantitative spectroscopy has been used to measure accurately the Doppler-broadening of atomic 
transitions in 85 Rb vapor. By using a conventional platinum resistance thermometer and the Doppler 
thermometry technique, we were able to determine Jzb with a relative uncertainty of 4.1 x 10 -4 , and 
with a deviation of 2.7 x 10 -4 from the expected value. Our experiment, using an effusive vapour, 
departs significantly from other Doppler-broadened thermometry (DBT) techniques, which rely on 
weakly absorbing molecules in a diffusive regime. In these circumstances, very different systematic 
effects such as magnetic sensitivity and optical pumping are dominant. Using the model developed 
recently by Stace and Luiten, we estimate the perturbation due to optical pumping of the measured 
ks value was less than 4 x 10 -6 . The effects of optical pumping on atomic and molecular DBT 
experiments is mapped over a wide range of beam size and saturation intensity, indicating possible 
avenues for improvement. We also compare the line-broadening mechanisms, windows of operation 
and detection limits of some recent DBT experiments. 



I. INTRODUCTION 

An ability to measure absolute thermodynamic tem- 
perature is currently confined to an exclusive list of na- 
tional standards laboratories. Such measurements re- 
quire primary thermometry expertise and highly spe- 
cialised equipment, making it impractical for wider-scale 
use. Instead, for practical purposes, temperature scales 
such as ITS-90 [l|, are used as an approximation to the 
true thermodynamic temperature, even though they are 
known to have relative uncertainties at the 10 -4 level Q. 
Accordingly, there is now a global push to develop more 
convenient methods of primary thermometry that could 
make thermodynamic temperature measurements more 
broadly available [3] . A second motivating factor driving 
recent renewed interest in primary thermometry is the 
call by the Bureau des Internationale de Poids et Mea- 
sures' to remeasure the Boltzmann constant, &b, using a 
wide range of techniques in preparation for the redefini- 
tion of the kelvin in 2011 0, 

At present, the recommended value for ks is derived 
primarily from one acoustic measurement performed by 
Moldover et. al. [6]. It would be preferable if fc# could 
be remeasured using different methods with compara- 
ble uncertainties so that any systematic errors can be 
identified. Recently, Borde suggested that the known Q 
temperature dependence of the spectral linewidth of ab- 
sorption features, previously exploited for measurments 
of stellar atmospheric temperatures [H, can be used as 
a new technique for accurate primary thermometry [9|. 
This Doppler-broadening thermometry (DBT) method 
has since been experimentally realized and while cur- 
rently less precise than other primary methods such as 
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acoustic or dielectric constant gas thermometry, it suf- 
fers from very different types of systematic errors [Toj| . In 
contrast to those techniques which measure the macro- 
scopic quantity RT, the DBT method can directly probe 
the microscopic thermal energy ksT by measuring the 
characteristic spectral width of molecular /atomic transi- 
tions, thereby avoiding the uncertainties in the Avogadro 
constant Na- 

The first DBT measurements performed by Daussy et. 
al. and Casa et. al. have used molecular gases at 1-130 
Pa and determined ks with a relative uncertainty 
[l3[ of order 10 -4 . More recent work has improved the 
statistical uncertainty [l4| to 4 x 10" 5 . In this paper, 
we present results from an atomic Rubidium (Rb) DBT 
system at 3 x 10 -5 Pa. The use of an alkali metal atom 
absorber presents some distinct spectroscopic differences 
that allow us to explore new experimental regimes within 
the DBT method. In doing so, we have identified some 
advantages in using an atomic system and encountered 
new challenges whose resolutions demand a deep under- 
standing and detailed unification of light-matter interac- 
tions and gas dynamical theories [l5j|. Whilst this work 
has been motivated by primary thermometry, we believe 
that it has wider implications in all fields requiring high- 
resolution or precision, quantitative spectroscopy. 

The primary difference between a low-pressure atomic 
vapor system, like Rb, and molecular experiments is that 
atomic motion is effusive, so collisions are extremely rare. 
An advantage of moving to this effusive regime is the 
avoidance of pressure-induced systematic changes, such 
as collisional line-shape perturbations [16| . This removes 
the need to extrapolate results to an equivalent zero- 
pressure value. Further, by using a sealed vapor-pressure 
reference gas cell, we avoid any unintended pressure drifts 
associated with a more sophisticated variable-pressure 
gas chamber. In this study, we chose Rb which has a 
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high optical cross section, ensuring that the absorption 
signal-to-noise is satisfactory in spite of the very low va- 
por density. This choice also gives the ability to operate 
in a range where high quality lasers are readily available 
and where silicon detectors, with their extremely good 
linearity [17j and quantum efficiency, can be used. At 
room temperature, it was possible to use a smaller cell 
than those in the molecular experiments for the same 
absorption depth, which was convenient for thermal con- 
trol. 

On the other hand, there are also disadvantages to us- 
ing atomic vapors. Care is required to ensure low levels 
of residual magnetic fields which can lead to the broad- 
ening of spectral lines by lifting the degeneracy of Zee- 
man levels. In the transition used in this experiment 
(D2 line of 85 Rb), optical pumping between unresolved 
hyperfine transitions contributed the greatest potential 
for systematic error by perturbing the lineshape. In ad- 
dition, the comparability of the natural lifetime, beam 
transit time and the Rabi period in our situation leads 
to optical pumping effects that requires either an accu- 
rate model of the gas kinetics [HI or the use of a low 
power probe beam to avoid lineshape perturbations. We 
followed the second route here and show here that the re- 
sulting perturbation is at the 4 parts-per-million (ppm) 
level. We note, however, that it should be possible to op- 
erate at higher power levels and use the model of Stace 
and Luiten [l5[ to correct the results. 

The remainder of this paper is organised as follows. 
In Sec. II, we describe the experimental setup and data 
acquisition for a Rb DBT experiment. The theoretical 
absorption coefficient lineshape for unperturbed atoms 
and the derivation of a suitable lineshape profile for ex- 
tracting the quantity ksT from the acquired spectra are 
then described in Sec. III. Our experimental results are 
presented in Sec. IV, followed by a discussion about sys- 
tematic uncertainties in Sec. V. In Sec. VI, we discuss 
the relative merits and challenges of the atomic approach 
in more detail. We will also discuss how the use of atomic 
vapors might be extended to a wider temperature range. 



II. METHOD AND EXPERIMENT 

A schematic of the experimental setup used in this 
work is shown in Fig. [TJ The ksT product was de- 
termined from spectroscopic measurements of unresolved 
hyperfine transitions (from the 5 2 Si/2, F = 2 — >• 5 2 P3/2, 
F = 1,2,3 states) in the D2 line of 85 Rb at 780 nm. 
A custom-built extended cavity diode laser (ECDL) was 
locked to a tunable optical Fabry-Perot cavit y (O C) us- 
ing the Pound-Drever-Hall (PDH) technique [18|, which 
suppressed acoustic and other fluctuations of the laser 
frequency. We tuned the laser's wavelength by chang- 
ing the length of the cavity. A sweep generator (Sweep) 
was used to displace one of the OC cavity mirrors by 
driving a piezo-electric stack actuator (PZT). A step- 
wise sweep was used to reduce the asynchronicity be- 




FIG. 1. (Color Online) A schematic diagram of the optical cir- 
cuit and electronic control systems. ECDL: extended cavity 
diode laser; LCD: variable optical depth liquid crystal dis- 
play; OC: optical cavity; PZT: annular piezo-electric stack; 
PDH: Pound-Drever-Hall top-of-resonance feedback control 
loop; RL: reference laser; VNDF: variable neutral density fil- 
ter; BE: beam expander lenses; PRT: ITS-90 calibrated plat- 
inum resistance thermometer; TI: thermal isolator; (F)PD: 
(fast) photodetector. 



tween measurements of the temperature, frequency and 
absorbed power. An LCD variable optical attenuator was 
used to keep the intensity transmitted through the cav- 
ity, which was monitored by the photodetector PD1, at 
a constant. This suppressed variations in the amount of 
optical pumping associated with power changes occurring 
throughout the scan. 

The optical frequency difference between the probe and 
a reference laser (RL) was measured using a fast pho- 
todetector (FPD) and a high bandwidth counter. The 
reference beam was derived from a tunable Ti: Sapphire 
laser locked to a temperature-controlled, ultra-low ex- 
pansion optical cavity. The stability of the reference laser 
(~1 kHz) contributed negligible uncertainty to the opti- 
cal frequency measurement. A variable neutral density 
filter (VNDF) allowed measurements to be made at inci- 
dent probe powers below that set by the LCD intensity 
control system, without affecting intensity or frequency 
stabilities. A photodetector (PD2) located just before 
the absorption cell was used to correct for any residual 
probe power changes. A beam expander (BE) was used 
to give the largest beam possible through the cell (2 cm 
in diameter) to maximise the input power for a given in- 
tensity to achieve a high signal-to- noise ratio. The power 
transmitted through a 10 cm long Rb absorption cell was 
measured on a third photodetector (PD3). The Rb cell 
was encased in a passive thermal isolator (TI) whose tem- 
perature was monitored using a using a platinum resis- 
tance thermometer (PRT) calibrated to ITS-90 with 30 
mK uncertainty. The maximum temperature difference 
across the isolator was measured to be less than 7 mK 
using auxiliary thermometers (not shown). The thermal 
relaxation time constant of the isolator (~2 hours) was 
much larger than the time required for a single sweep (10 
mins). 
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FIG. 2. (Color online) The lower panel shows a typical spec- 
trum (blue circles) obtained from Rb vapour at 295 K (~550 
points) and the fit (red) using Eqn. (|2]). The Lorentzians 
(green) indicate both the location of the underlying hyperfine 
transitions and their relative strengths (in arbitrary units). 
The upper panel displays the residuals to the fit. 



The outputs of PD2 and PD3 that monitored the in- 
cident probe and the transmitted powers, respectively, 
were recorded using digital multimeters with 5-digit res- 
olution. A third multimeter was used to make four-wire 
resistance measurements of the PRT. A computer-based 
data acquisition system centrally controlled and logged 
the three digital multimeters and the high bandwidth 
counter outputs. Before each set of measurements at a 
fixed probe power level, the DC voltages on PD2 and PD3 
caused by stray lights, were recorded in the absence of 
the probe beam. These backgrounds were removed from 
the detector readings before reconstructing the Doppler- 
broadened absorption spectrum by dividing the output 
of PD3 by that of PD2. 



III. EXTRACTING k B T FROM SPECTRA 

In a low density gas cell, the absorption coefficient of 
an isolated spectral line is the convolution between a 
Lorentzian and a Gaussian function, known as a Voigt 
profile [7]. The Lorentzian component can be written as 
L{y) = (1 + [{y — ^o)/r] 2 ) -1 , where v is the optical fre- 
quency and uq is the atomic transition frequency. This 
component arises from the finite lifetime of the upper 
state and has a full-width-at-half-maximum (FWHM) 
of T = 1/(27tt), where r is the total upper state life- 



time. Various phenomena including collisional, transit 
time and power broadening can perturb the width of this 
Lorentzian. It was important in this experiment to keep 
these perturbations at negligible levels. 

The Gaussian component has the form of G(y) = 
exp(— [(y — vq)/ Ai^d] 2 ), where the 1/e-half- width of the 
Gaussian component is related to the thermal energy 
ksT as: 



k B T 



_ mc 2 f Az/p V 



(i) 



and arises from the thermal motion of the atoms. Here 
c is the speed of light, m is the mass of the absorb- 
ing atom and v$ is the absolute transition frequency. 
For Rb, the latter two parameters are known with rel- 
ative uncertainties of 5 x 10 -8 and 5 x 10 -11 , respec- 
tively pj| [20|, which are negligible for this experiment. 
Whilst it was possible to compute the absorption coeffi- 
cient V(y, Avjj) = L(y) ®G(y) by performing a convolu- 
tion integral, this was a computationally inefficient proce- 
dure. Instead, a numerical routine described by Humlicek 
[2lj implemented in the scientific data analysis program 
Igor Pro [22I] was used to generate the Voigt function. 
The relative accuracy of this algorithm was better than 
3 x 10 -5 and no consideration of these numerical errors 
was required in the following analysis. 

The observed transmission spectrum was determined 
by Beer's Law [7], T(y) = exp(— V(v, Az/d)L), where L 
is the optical path length of the probe inside the ab- 
sorbing gas. Since this experiment concerns frequency 
differences, it was convenient to shift the frequency ori- 
gin by subtracting away the absolute optical frequency 
of the 87 Rb F = 2 — >• 1 hyperfine transition. We denote 
these frequency differences with / to avoid confusion with 
the absolute frequencies v. We performed least-squares 
regression of the transmission spectra to Eqn. (|2]) to de- 
termine Az/p. 



T(f)=AfyBexp[-Cj2SiVi(f-fi 



-/cav,Az/ D )j(2) 



Here, T(f) is the sum of three, unresolved hyperfine tran- 
sitions, indicated schematically in Fig. [2j the Vi(f — 
fi — /cav) are the Voigt functions corresponding to each 
hyperfine transition, each centered on frequency fi and 
with relative strength Si m , and A, B and C are three of 
the five parameters adjusted in the regression. The Si 
were fixed at values derived from the Clebsch-Gordan 
coefficients while the U were fixed at values determined 
by Arimondo et. al [20[ Of the adjustable parameters, 
A and B account for, respectively, a residual background 
and imperfect background normalisation (to a value of 1) 
of the spectra. C is the on-resonance absorption depth 
determined by the number density and the cell length. 
The two other adjustable parameters, / cav and Avd were 
contained within, and were common to the three Voigt 
functions, with / cav allowing for the arbitrary (but fixed) 
frequency offset of the reference laser from the lowest 
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FIG. 3. (Color Online) The ( X 2 

change in Lorentzian and Gaussian widths, denoted F 
and Aisd — Az/D, op t respectively, show that there is a high 
degree of correlation between these fit parameters. The sub- 
scripts opt denote the optimum quantities determined from 
the least-squares fit, which were 2.75 MHz and 308.0 MHz, for 
the Lorentz and Gaussian widths, respectively. The thicker, 
red contour is generated from the covariance matrix and is 
intended to guide the eye. 



TABLE I. Values of the fixed constants used in Eqn. [2] for the 
relative transition strengths, Si, and their frequency offsets 
[20l ], fi. The relative transition strengths are exact numbers, 
whilst there is a relative uncertainty of order 1 x 10 -5 asso- 
ciated with the differences between the hyperfine transition 
frequencies. 
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frequency hyperfine line. The Boltzmann constant was 
obtained from the Gaussian width Az/p using Eqn. (pQ). 

In principle, it is possible to perform the regression al- 
lowing every parameter in Eqn. [2]to be adjustable. How- 
ever, in this case, the returned parameters from the least- 
squares fit show strong correlations between a number of 
the free parameters. This leads a larger uncertainty range 
in those parameters. In particular, we find strong corre- 
lations in the Si, and also in T and Aujj which have a 
correlation coefficient of -0.96. Operating in an intensity 
range where it is possible to fix the Lorentz width and 
the strengths Si to known values, we avoid an artificial 
inflation of the uncertainties associated with the fitted 
Doppler width due to these correlations. Our theoretical 
model (see Section |V|) guides us on the upper intensity. 



IV. RESULTS 

The results from several ks determinations are shown 
in Fig. Has a function of incident probe power. These 
data were produced from 24 spectra taken at intensi- 
ties at least 500-fold below the saturation intensity [23| 
(^sat ^20 W/m 2 ). The total uncertainties at each power 
level are indicated by the error bars and were determined 
from the standard deviation of the data sets taken at the 
same power level. Our final determination of the Boltz- 



FIG. 4. (Color Online) The value of Jzb (dots) extracted 
from fitting spectra measured at several probe powers. The 
weighted mean (blue dashed) and weighted one-standard- 
deviation band (blue solid) are shown. The green, dot-dashed 
line is the CODATA value. 



mann constant, based on a weighted mean of the mea- 
surements at the various probe powers, is 

k B = 1.38104(59) x 10~ 23 J/K. 

The uncertainty in this measurement was limited pri- 
marily by residual amplitude noise at frequencies higher 
than the bandwidth of the LCD feedback control loop. 
This white noise was exacerbated by a lower than ex- 
pected common-mode rejection of amplitude noise re- 
sulting from some asynchronicity in the data acquisition 
scheme. 

The broadening of the Lorentzian linewidth by mag- 
netic fields contributed a smaller systematic uncertainty 
of 9.8 x 10 -5 , which was summed in quadrature give a 
total relative uncertainty of 4.1 x 10 -4 . All other sys- 
tematic uncertainties were kept at negligible levels and 
will be discussed in the following section. However, as 
the random uncertainties and measurement timing are 
addressed in future work, these will become increasingly 
important. Our value of fc# has a relative deviation of 
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2.7 x 10" 4 from the current CODATA value. To our 
knowledge, this is the most accurate demonstration of 
primary thermometry using spectroscopy of atomic ab- 
sorbers. 



V. ESTIMATION OF SYSTEMATIC 
UNCERTAINTIES 

The greatest source of systematic uncertainty in this 
experiment was caused by the Earth's magnetic field. 
This external field lifted the frequency degeneracy of the 
Zeeman sub- levels in the Rb hyperfine transitions. The 
field along the cell was measured to be 0.21(1) G and nu- 
merical simulations shows that this causes an apparent 
broadening of the Lorentzian linewidth by 0.5%, which 
in turn, contributes a relative uncertainty of 9.8 x 10 -5 
in the measured value of ks- Whilst magnetic shields 
can be used to reduce the field strength by a factor of 
40 000 times, this was not attempted because the effect 
was smaller than the observed random uncertainties of 
the experiment. 

The Rb spectra were recorded (non-contiguously) over 
a one month period. Whilst no drift would be expected 
from the use of a permanently sealed Rb vapor-pressure 
reference cell with no buffer gas, checks for long-term sta- 
bility were performed by periodically re-acquiring spec- 
tra at previously investigated probe power levels. We 
detected no long term drift larger than the short-term 
reproducibility of the experiment. A long-term drift of 
the reference laser's wavelength only changes the fitted 
offset frequency fo but does not change the Rb spectrum 
and is therefore inconsequential. 

We also considered the systematic uncertainty caused 
by holding T and Si constant. Fig [3] demonstrates 
that the minimum of the y 2 surface lies along (T — 
r opt )/(A^£> — Az/D j0pt ) ~ 1, so the propagated relative 
uncertainty in Aujj due to the uncertainty of T is re- 
duced by a factor ~ Avjj/T « 100. Since the unper- 
turbed Lorentzian width is known [24| with a relative 
uncertainty of 3 x 10 -4 , this contributes a negligible un- 
certainty in Aud. Similarly, there is no uncertainty as- 
sociated with the transition strengths because they are 
exact numbers. However, these assumptions are valid 
only in the limit of zero probe beam power because the 
measurement process itself can perturb the actual val- 
ues of T and Si away from their theoretical values. For 
example, excess probe beam power will lead to power- 
broadening of the underlying Lorentzian linewidth, in- 
troducing an intensity dependence in the fitted Doppler 
width and, thus, in fc#. Similarly, any optical pumping 
between the hyperfine states that produces alignment or 
polarization of the sample will perturb the effective ratios 
Si/Sj. We did indeed observe such effects at intensities 
greater than ~1 W/m 2 (at least 30 times larger than in- 
tensities reported here), which is still more than an order 
of magnitude below the saturation intensity (23| . 

We have developed a new theory [15] to explain the 



TABLE II. Different experimental regimes accessed by vari- 
ous choices of absorbers. The relative noise reported in the 
bottom two rows are given for a 1000 s bandwidth. 





Rb 


NH 3 [11] 


CO2 [12] 


Doppler Width (MHz) 


308 


50 


275 


Pressure (Pa) 


10~ 5 


8 


130 


Pressure Broadening (MHz) 


0.02 


1 


4 


Optical Pumping, & 


0.021 


8xl0" 8 


3xl0" 7 


Probe Power (fiW) 
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50 


1000 s Shot Noise Floor (ppm) 


0.7 


0.8 


1.2 


1000 s Achieved Reproducibility (ppm) 


4000 


1000 


200 



cause of this early onset of optical saturation behaviour 
which allows us to confirm that a negligible amount of 
optical pumping occurred in the experiments reported 
here. In contrast to previous models [l7], US HH that 
approximate the atom as a two-level system under uni- 
form illumination, we developed a theory for multi-level 
atoms probed by a beam with finite spatial extent [l5| . 
This better captures the physical properties and situation 
of this Rb experiment. We find that the saturation-like 
behaviour comes from an optical pumping process that 
evacuated the population in the laser-coupled ground 
state, and not from a Rabi-flopping mechanism responsi- 
ble for saturation in a two- level system 0, H^-We, there- 
fore, quantify the perturbation away from thermal equi- 
librium by introducing a figure-of- merit, denoted de- 
fined as the relative depletion of the laser-coupled ground 
state, i.e. 

& = 1 - — (3) 
Pi,th 

where pi,th and pi are the thermal and pumped frac- 
tional occupation of the ground state. In the absence of 
optical pumping (pi = pi,th), & = 0; whilst conversely, 
when the ground state is fully depleted, (pi = 0) and 
& — 1. Using the full model of Stace and Luiten [lEj . 
we computed the figure-of- merit to be & ~ 0.021 for the 
largest powers in this experiment. In Appendix [Aj we 
show how & may be calculated approximately using a 
simple method based on conservation of energy consid- 
erations. At & = 0.021, the perturbation in ks due to 
optical pumping is 4 ppm (see Appendix [B]), which is 
well below the current statistical uncertainty. We also 
showed experimentally that we were sufficiently close to 
the zero-probe-power limit by demonstrating the absence 
of any statistically significant slope in the deduced value 
of ks as a function of power on Fig. |U 



VI. COMPARISON TO MOLECULAR DBT 

In this section, we will highlight the differences be- 
tween a low-pressure atomic approach and high-pressure 
molecular DBT experiments by examining reported re- 
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suits, simulated spectra and the results presented in Sec- 
tion IIV[ Table [IT] shows a comparison of some recent 
DBT experiments using molecules and atoms. 



A. Shot-noise limits 

Using the probe powers shown in Table HU row 5, 
we simulated shot-noise limited Doppler-broadened spec- 
tra consisting of 1000 absorbance samples spanning /q ± 
3Au]j. An approximate estimate of the shot-noise lim- 
ited relative uncertainty of ks was found by fitting these 
spectra to our model. The results are surprisingly similar 
for all experiments (row 6). The results depend mostly 
on the total scan time and not details such as the fre- 
quency spacings between points. Despite the high CO2 
probe power reported in Ref. [Hj, the shot-noise limit 
is similar to the other experiments because of the low 
absorption depth in that experiment. The last row of 
Table Ull shows that the approximate actual repeatability 
calculated from reported results, when scaled to 1000 s 
time scale, is significantly worse than the shot-noise limit 
in all cases. Djerroud et. al. [14] have reported recent 
improvements to the work of Daussy et. al. by achieving 
a reproducibility of 100 ppm at 1000 s. In our case, we 
found that background light levels have a significant im- 
pact on experimental repeatability, so we intend in future 
experiments to employ synchronous detection of a mod- 
ulated carrier. Preliminary results suggest this improves 
the repeatability to 150 ppm. 




10~ 6 10~ 5 10~ 4 10~ 3 10~ 2 10~ 1 

(n/if 

FIG. 5. (Color Online) A contour plot of & computed from 
the model of St ace and LuitenJT^], as a function of normal- 
ized probe intensity (Q/r) 2 and normalised coherent evolu- 
tion time (Tt). The red and blue dots indicate the condi- 
tions of this and the molecular DBT experiments, respec- 
tively, demonstrating the exploration of vastly different spec- 
troscopic regimes. The dot-dashed lines are contours of con- 
stant input power. The solid red line indicates the experi- 
mental conditions that would lead to a 1 ppm perturbation 
in ks- 



B. Optimization of Probe Beam Power and 
Geometry 

Current a ppr oaches in DBT experiments range in 
pressures [l2[ Tlij from 10 -5 to 100 Pa. In the Rb exper- 
iment described here, the mean free path is 120 m and 
the gas dynamics is governed by an effusive flow of atoms, 
i.e. the system is essentially collisionless, leading to neg- 
ligible pressure-induced self-broadening Q. In contrast, 
the molecular DBT experiments are dominated by colli- 
sions and is in the diffusive regime, causing ~ 0.01 Aujj of 
pressure-related broadening [ll|,[28|. Moreover, the phys- 
ical details of the molecule-molecule collisions are critical 
in determining the observed lineshape profile. For exam- 
ple, the degree of elasticity of molecular collisions per- 
turbs the lineshape away from a Voigt profile [l6j|. This 
requires higher-order forms such as the Galatry profile 
and a greater number of adjustable parameters required 
for the regression of spectra. 

In contrast to the collisional effects of the molecu- 
lar DBT systems the critical linewidth perturbation for 
vapours with strong light interaction is associated with 
optical pumping. This is clearly demonstrated in row 4 
of Table [III To illustrate the very different operational 
spectroscopic re gim es interrogated in this work and by 
Daussy et. al. [Ill, H and Casa et al BO, Fig. 



[5] shows a contour plot of & as a function of the non- 
dimensional quantities Ft and (Q/T) 2 . The vertical axis 
corresponds to the ratio of the coherent interaction time, 
t, to its lifetime (1/T): for an effusive vapor like Rb, t is 
the average beam transit time while for a diffusive gas, 
like those used by the molecular experiments, t is the 
mean time between collisions. The horizontal axis is pro- 
portional to the beam intensity, and is expressed in terms 
of the Rabi frequency, ^, and natural decay rate, T. The 
conventional intensity can be related to these parame- 
ters by the expression / = 2(£7/T) 2 / sat , where J sat is the 
conventional 2-level saturation intensity [23]. The & con- 
tours, which were calculated with the model of Stace and 
Luiten[15|, can be closely approximated by the simpler 
expressions: 

^ f (nt/2) 2 im»2 

~ \ n 2 t/(2T) if Tt < 2 " V J 

Thus, contours of constant & have a slope of — 1 
decade/decade and —1/2 decade/decade when Tt ^> 2 
and Tt <C 2, respectively. Eqn. 2] was derived using a 
model of a three level atom with two ground states and 
an excited state. The scales of Fig. [5] obscure the details 
of the more complex case Tt « 2 , which are unimportant 
to this wide range overview. 
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In Fig. [5j the solid red line indicates the parameters 
where there would be a 1 ppm perturbation to the deter- 
mination of ks (derived from the analysis in Appendix 
|B|) . For collisionless gases, t is proportional to the beam 
radius r and (Q/T) 2 oc r -2 . Therefore, changing the 
probe beam diameter for a fixed input power lead to lines 
with a slope of —1/2 decade/decade. Two such contours 
(dot-dashed) are drawn through the location of the Rb 
experiment reported here as well as for the molecular ex- 
periments. Contours of constant absorbed power follow 
lines of constant & since the absorption coefficient is ap- 
proximately proportional to the population perturbation 
in the ground state, i.e. 



V(f u Au D ) oc Pl =pi,th(l -^). 



(5) 



It is clear from Fig. [5] that there are two regimes 
of coherent interaction time, probe intensity and opti- 
cal pumping which are delineated by the line Ft = 2. 
For Tt <C 2, the contours of & and those for fixed in- 
put power are parallel, demonstrating that the measured 
results will be independent of the probe beam geometry. 
However, this does not hold when Tt ^> 2, indicating 
that the result will depend on probe beam shape. In this 
situation, Fig. [5] guides the experimenter into an opti- 
misation of probe beam radius and intensity for a given 
perturbation to fc#. By increasing the beam radius, it 
is possible to increase the total input power whilst re- 
maining on the red line. This is advantageous as it will 
lead to an improved SNR of the measurement for a given 
linewidth perturbation. The limit to this procedure is 
set by the size of the gas cell. For molecular experi- 
ments where the coherent interaction time is limited by 
the mean time between collisions (and not the beam tran- 
sit time), an increase in probe power will require a cor- 
responding increase in pressure to preserve the amount 
of optical pumping. Of course, this compensation may 
be forgone given the much smaller magnitude of optical 
pumping in this regime. However, the small optical cross 
sections of the molecular transitions which reduce the im- 
pact of optical pumping effects give rise to other issues 
associated with SNR and pressure-broadening. 



C. Absorption Depth and SNR 

Whilst Fig. [5] demonstrates the impact of experimen- 
tally controllable variables such as the coherent interac- 
tion time (which can be adjusted by changing the beam 
radius in an effusive gas system, and by the pressure in 
a diffusive gas) and the probe beam intensity, there are 
other considerations pertaining to the atomic/molecular 
properties that independently affect SNR. One such con- 
sideration is the optical depth. Numerical simulations 
show that in the presence of white amplitude noise (such 
as detector noise) there is an optimum optical depth, 
aL « 3 at which the variance in the Doppler width ob- 
tained from repeated simulations is minimised. In terms 




200 300 400 500 600 700 800 9001000 

Temperature (K) 

FIG. 6. (Color online) Approximate temperature windows of 
operation for various vapor- pressure cells, of length restricted 
to 10 cm. The species in each cell is shown on the vertical 
axis. The temperatures of selected ITS-90 fixed points are 
indicated by the solid vertical red lines. The subscripts TP, 
MP and FP denote Triple Point, Melting Point and Freezing 
Point, respectively. 



of fundamental parameters, the optical depth can be ex- 
pressed asp7| 



aL = cfqPqL oc TX PL, 



(6) 



where <to, po> P and L are the optical cross-section, num- 
ber density, pressure of the absorber and optical path 
length, respectively. We can more clearly see the trade- 
offs between the parameters that affect the optical depth 
by writing an "equation of state" for the optical depth us- 
ing the conditions of the Rb experiment as the reference 
values, 



aL 
07 
(A/nm) 2 
780 2 



T/s" 1 \ / P/Pa 
4 x 10 7 ) V 3.5 x 10- 5 



L/m 

IF 1 



(7) 



To achieve comparable SNR to atomic absorbers in which 
T is large, molecular-based approaches must compensate 
by increasing the product of P and L because of the much 
longer lifetimes. Increasing either of these parameters 
can lead to the potential for unwanted inaccuracies in the 
experiment relating to pressure broadening or thermal 
gradients in the apparatus. 



D. Windows of Operation 

In this experiment, we used a commercially available 
Rb vapor cell of a standard length (10 cm), which gave 
an absorption depth of ~ 50% at room temperature. 
However, operating along the metal's sublimation curve 
means that the equilibrium vapor density increases ex- 
ponentially with temperature, and thus the absorption 
depth is much more sensitive to temperature for atomic 
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systems than in the molecular experiments. The range of 
temperatures over which we can use Rb as a thermomet- 
ric substance is comparatively smaller than for molecules. 
The lower temperature bound is provided by the dimin- 
ishing absorption depth resulting from the exponential 
decrease in vapor pressure. Conversely, the upper tem- 
perature bound is set when the vapor density becomes so 
large that the probe beam is completely absorbed over 
a small range of frequency detuning. Fig. [6] shows ap- 
proximate temperature windows over which a selection 
of metal gases at their vapor pressure might be able to 
operate [23|, l30l - [36| . At the upper and lower temperature 
bounds of each window, the irreproducibility of the fitted 
Doppler width in repeated simulations is double that for 
the optimal optical depth (i.e. aL « 3). This condition 
restricts the optical depth to a range 0.5 < aL < 90. We 
have restricted attention to 10 cm long cells for conve- 
nient thermal control. Some fixed points of the ITS-90 
[lOj scale are also shown, demonstrating that metallic 
vapor-pressure DBT can be used to verify the thermody- 
namic temperature at those fixed points. 

For a particular choice of absorber, it is also possible to 
fill the reference cell using an external atomic reservoir at 
a lower temperature, resulting in a lower vapor number 
density. Once the vapor cell is sealed, no further increase 
in number density is possible, thereby restricting optical 
depth to an optimum level when at high temperatures. 
This would permit the exploration of temperatures and 
pressures away from the sublimation phase boundary. 



VII. CONCLUSIONS 

Doppler broadening thermometry provides a signifi- 
cantly different approach to the determination of the 
Boltzmann constant in preparation for the re-definition 
of the kelvin. Whilst previous DBT experiments have 
studied molecular absorbers at up to 10 Pa in a diffusive 
regime, we have used Rb at 3 x 10 -5 Pa where the gas dy- 
namics are effusive. Our approach avoids problems such 
as pressure broadening and poor signal-to- noise. Fur- 
thermore, the Rb cell is compact in size, enabling more 
convenient temperature control. We used this system to 
determine ks with a relative uncertainty of 4 x 10 -4 . 
The present experiment was limited by amplitude noise 
in the probe beam, which in future work will be overcome 
by using a control loop with larger bandwidth. We esti- 
mate that all current DBT experiments are far from the 
shot noise limit and that the use of synchronous detection 
methods could address this situation. 

We also compared sources of systematic uncertainty 
for DBT experiments in the diffusive regime of molec- 
ular DBT and the effusive regime of atomic DBT. For 
atomic vapors such as Rb, the shift in the determined 
value of ks due to the Earth's magnetic field is approxi- 
mately 100 ppm. Optical pumping effects are determined 
by the combination of beam transit time, upper-state 
lifetime and probe intensity. In this work, the equilib- 



rium ground-state population was perturbed by 0.2% 
which shifted the measured ks value by 4 ppm. The 
small cross-sections of the transitions used in molecular 
DBT experiments means that the optical pumping effects 
are relatively small (perturbation to equilibrium ground- 
state population of order 10 _T ). However, these small 
cross-sections mean that an equivalent signal-to-noise in 
molecular experiments requires long path lengths and/or 
higher pressures, both of which have associated system- 
atic effects on measurements of ks- The method devel- 
oped here for quantifying optical pumping effects can be 
used to search for candidate absorbers optimally suited 
to DBT experiments in either the diffusive or effusive 
regimes. 

Appendix A: Approximate method for calculating & 

Whilst the approach of Stace and Luiten [l5| cap- 
tures the details of the gas-dynamical interaction between 
atoms and the beam volume, which can be used to esti- 
mate the population perturbation using no phenomeno- 
logical parameters, it relies on a knowledge of the dipole 
matrix element in order to calculate the Rabi frequency. 
It is useful to compare this against a simpler, independent 
analytic calculation. We provide such a calculation here, 
based on energy conservation and using an a posteriori 
knowledge of the optical depth. 

Assuming that the input intensity is sufficiently weak 
so that an atom never absorbs more than one photon 
during a single transit of the beam, the rate of atoms 
scattered out of the laser-coupled ground state is approx- 
imately the same as the rate of scattered photons. This 
is given by 

R = a^-(1-I/I ), (Al) 

where a = irr 2 is the cross-sectional area of a beam with 
radius r, is the photon energy and / and Iq are the 
on-resonance output and input intensities. 

Using kinetic theory [37j, the flux of atoms crossing 
the beam surface (of area A = 2nrl, where I is the cell 
length) is 

$A=\vp Q -^—A, (A2) 

where v = ^2kT/ m is the mean atomic speed and po is 
the total number density. A factor of T/Avd is used to 
count only those atoms which are on-resonance. 

The value of & (Eqn. [3] from Section [Vj) is then the 
ratio of the rate of scattering to the rate of new atoms 
impinging on the beam, i.e. 

& = R/{?bA) 

« 0.022 for this experiment. 
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FIG. 7. (Color Online) An example of the unperturbed ab- 
sorption coefficient (ao(A)) and the additive perturbation 
term (r(A)), plotted for v « 30MHz/3MHz = 100. The per- 
turbation term is vertically expanded by 10 5 for clarity. Note 
that v in this Appendix has a different meaning to that used 
in the preceding sections of this paper. 



This approximation agrees with the more detailed anal- 
ysis, differing only by 10 -3 . 



Appendix B: Linewidth Perturbation 

Using the approach and notation of Stace and Luiten 
[Hj , the general expression for the absorption coefficient 
(Eqn. (22) of Ref. |l5j) can be expanded (Appendix Al) 
in the case of small optical pumping. The expression for 
a perturbed absorption coefficient lineshape is 



a 




KcP\Z 

2uo Pi 8v 

P£ -P{' 4(A/v) 2 -3 



1 



M 

A 2 



V(A,u) 



Pi 64^3 
a (A) + er(A), 



■ exp 



1/4 



(Bl) 



where ceo (A) is the unperturbed lineshape, r(A) is de- 
fined as the exponential term in the parentheses and A 
is the detuning in units of the Gaussian width. Fig. (0 
shows a plot of the unperturbed lineshape along with 



the perturbation term r(A). We performed an expan- 
sion in the small parameter e = (P^ — P")/Pi, where 
pn = d 2 p i /dA 2 the derivative is with respect to the de- 
tuning A. The meanings of the remaining variables in 
this section are as defined in Ref. [l5|. Written in this 
form, it is clear that the absorption coefficient in the pres- 
ence of some optical pumping is a Voigt function, denoted 
V(A, v) with an additional perturbation term. 

Consider the expansion of a characteristic width for 
the perturbed absorption coefficient, 



Ai/ e ,pert( e ) ~ A 1 / e?pert (0) + 



dA 1/( 



e,pert 



5e, (B2) 



/ejpertvw ~ ^l/e,pertv^ / 1 j 

such that «(A 1 / e pert )/a(0) = 1/e. Differentiating this 
definition and using A 1 / e pert (0) = Ai/ e , we obtain the 
perturbation to the linewidth 



dA 1 



/e,pert 



de 



e=0 



a (0)[r(0)/e-r(A 1/e )] 
[a (0) + er(0)]K(A 1/e ) + er'(A 1/e )] 
a (0)[r(0)/e-r(A 1/e )] 
a a / (A 1/e ) 
1/e 1 



16^ 3 a 7 (A 1/e ) 
1 



■2 " 



(B3) 



To evaluate Se, we have used the low probe power ap- 
proximation of Stace and Luiten (P3 w << Pi) and 
A 1/e « v = Av D /T w 100 » 1. We also used the 
approximate form for the steady-state population den- 
sity Pi that ignores atomic coherence. This quantity is 
maximal on resonance (i.e. when A = 0): 



Max[P{'] = P{'(0) 



d 2 



1 f}* 2 r 

-(l-exp(- i - aa ^ )) 



< 



dA 2 \2 K ^ v 1 + 4A* 2 
-4**ir 2 exp(-**fi* 2 ) 
-4t*^* 2 

-0.4 for this experiment. 



A=0 



Substitution of the relevant quantities into Eqn. (|B2|) 
gives 2 ppm perturbation to the linewidth and, therefore, 
4 ppm in the determination of ks for this experiment. 
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